Naturally light invisible axion and local Z i3 ® Z ?J symmetries 
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We show that by imposing local Z13 (g> Z3 symmetries in an SU(2) <g> U(l) electroweak model we 
can implement an invisible axion in such a way that (i) the Peccei-Quinn symmetry is an automatic 
symmetry of the classical Lagrangian; and (ii) the axion is protected from semi classical gravitational 
effects. In order to be able to implement such a large discrete symmetry, and at the same time allow 
a general mixing in each charge sector, we introduce right-handed neutrinos and enlarge the scalar 
sector of the model. The domain wall problem is briefly considered. 
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It has been known for a long time that the complex 
nature of the QCD vacuum solves the old U(\)a prob- 
lem [1,2] but also implies the existence of the OFF term 
in the classical Lagrangian [3] . This CP-nonconserving 9 
angle (which has also an electroweak contribution) should 
be smaller than 1CP 9 in order to be consistent with the 
measurement of the electric dipole moment of the neu- 
tron [4]. Why is CP only weakly violated in QCD? This 
is the strong CP problem and an elegant way to solve 
it is by introducing a chiral symmetry £/(1)pq [5]. This 
Peccei-Quinn (PQ) symmetry implies also the existence 
of a pseudo Goldstone boson — the axion [6]. The invisi- 
ble axion is almost a singlet (here denoted by (f>) under 
the SU(2) ® U(l) gauge symmetry [7]. In addition, since 
this particle can occur in the early universe in the form of 
a Bose condensate which never comes into thermal equi- 
librium, it can be a significant component of the cold 
dark matter if its mass is of the order of 10 -5 eV [8,9], 
which makes it important in its own right. However such 
a small mass for the axion can be spoiled by semiclas- 
sical gravity effects since it could induce renormalizable 
or non-renormalizable effective interactions which break 
explicitly any global symmetry, and in particular the PQ 
symmetry. The most dangerous of these effective opera- 
tors are of the form (f> D , where D is the mass dimension 
of the effective operator. In this case the axion could gain 
a mass which is greater than the mass coming from in- 
stanton effects and also the angle is not less than 1CP 9 
in a natural way. This can be avoided if the dimension of 
the effective operators is D ~ 12 [10,11]. Hence, unless 
D is high enough, invisible axion models do not solve the 
strong CP problem in a natural way. Thus the invisi- 



ble (very light) axion must be protected against gravity 
effects. 

Some years ago Krauss and Wilczek [12] pointed out 
that local symmetries can masquerade as discrete global 
symmetries to an observer equipped with only low-energy 
probes. We can see this as follows. Consider a U(l) 
gauge theory with two complex scalar fields r\ and £ car- 
rying charge iVe and e, respectively. It is possible, for 
instance, that r\ undergoes condensation at some very 
high energy scale v, but £ produces quanta of relatively 
small mass and does not condense. Before the conden- 
sation the theory is invariant under the transformation 
rj'[x) — exp [iNeA(x)]r)(x), = exp [ieA(x)]£(x) and 
A'^ix) = A ll (x) — d ll A(x), here is the U(l) gauge field. 
The condensate characterized by the vacuum expectation 
value (r)(x)) = v in the homogeneous ground state is in- 
variant only when A is an integer multiple of 2ir/Ne. 
This residual transformation still acts nontrivially on £ 
but now as a discrete Zn group. The effective theory well 
below v will be the theory of the single complex scalar 
field £, neither the gauge field nor the scalar rj will appear 
in the effective theory since these fields are very heavy. 
The only implication of the original gauge symmetry for 
the low energy effective theory is the absence of interac- 
tion terms forbidden by the Zn symmetry. For instance 
if there were more charged scalar fields in the theory, 
the discrete symmetry would forbid many couplings that 
were otherwise possible. Hence, since no processes like 
black-hole evaporation or wormhole tunneling can vio- 
late a discrete gauge symmetry it means that local Z^ 
symmetries violate the no- hair theorem [13]. These sym- 
metries may be the consequences of gauge U(l)'s sym- 
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metries which are embedded in a larger grand unified 
theory or from a fifth dimension as it was shown for the 
case of U(l)y in Rcf. [14] or even in superstring inspired 
models [15-17]. 

In fact, it has long been known that superstring theo- 
ries exhibit PQ-like symmetries. It happens that in the 
reduction of a ten-dimensional theory to four dimensions 
an arbitrary large number of massless (up to the effect 
of the coupling FF) scalar four-dimensional fields, which 
cannot be gauged away, arise [18,19]. On the other hand, 
in superstring theory, also naturally local discrete sym- 
metries arise [15-18]. These facts are welcome since su- 
perstring theory is considered a good candidate for the 
unification of the four interactions. Hence it is reason- 
able to require that any global or local symmetry that 
is introduced in low energy models can be derived from 
the superstring theory. However, since we already do not 
know the details of the theory that describes quantum 
gravity, we will apply the ideas above in a more conser- 
vative scenario: a multi-Higgs extension of the standard 
model. 

We show below that in this situation, the axion, as a 
solution to the strong CP problem, can be achieved by 
introducing appropriated discrete symmetries and also 
the PQ symmetry arises as an automatic or accidental 
symmetry of the classical Lagrangian in the sense that 
it is not imposed on the Lagrangian but it is just a con- 
sequence of the particle content, renormalizability, and 
the gauge and Lorentz invariance of the model [10]. The 
automatic PQ symmetry has been also considered up to 
now only in the context of grand unified or supersymmct- 
ric theories [20] or, recently, in 3-3-1 models [21,22]. The 
same can be said with respect to the axion protection 
from gravity effects [11,16,17,22]. 

The model we will put forward can be viewed 
way of stabilizing the Dine-Fischler-Srednicki (DFS) ax- 
ion [7] . As is well known the DFS axion suffers from two 
difficulties: (i) the PQ symmetry has to be imposed by 
hand and (ii) the axion is not stable from the gravita- 
tional effects as discussed before. With respect to the 
first problem we expect that the PQ symmetry is not an 
ad hoc symmetry but an accidental one. It could be the 
consequence of discrete symmetries Zn which in addi- 
tion, if N is large enough, it can stabilize the axion from 
the semi-classical gravity effect. In this vain a large Zm 
symmetry solves both difficulties of the DFS axion model. 
The lowest order effective operator contributing to these 
parameters have D = 13 since operators with dimension 
up to D = 12 are forbidden by the local discrete sym- 
metry .Z13 £g> Z-j,- Moreover, after imposing these discrete 
symmetries the PQ chiral symmetry is automatically im- 
plemented in a model with SU(3) C <8> SU(2) L <g> U(1) Y 
gauge symmetry (a similar analysis for avoiding B- and 
L-violating operators in supersymmetric versions of this 
model was done in Ref. [23]). 

The axion is made invisible if it is almost a singlet 
under the gauge group and for this reason we add a sin- 



glet </> which gets a vacuum expectation value denoted by 
«0 [7] . The suppression of the explicit breakdown of the 
PQ symmetry by gravity is obtained by adding a large 
Zm symmetry among the matter multiplets, in such a 

way that effective operators like 4 }N ~ l /M^ ^ 4 , where 
Mpi is the Planck mass, are automatically suppressed. 
At the same time this local discrete symmetry makes 
the PQ symmetry an automatic symmetry of the clas- 
sical Lagrangian. For instance, a Z\% symmetry implies 
that the first nonforbidden operator is of dimension thir- 
teen and the first contribution to the axion mass square 
is proportional to w"/ M pi ~ 10~ 21 eV 2 or 10" n m 2 , if 
m a <~ Aq CD /u0 ps 10~ 5 eV is the instanton induced mass 
(we have used M Pl = 10 19 GeV and w = 10 12 GeV). 
We see also that the naturalness of the PQ solution to 
the ^-strong problem is not spoiled since in this case we 
have 9 c fi cx /Mj?J _4 Aq CD [24] and it means also that 
6 cS oc 10- 11 for N = 13. 

In order to be able to implement this rather large 
symmetry the representation content of the model is 
augmented: we add right-handed neutrinos and several 
scalar Higgs multiplets. Presently, right-handed neutri- 
nos may be a necessary ingredient of any electroweak 
model if, as indicated by recent experimental solar [25] 
and atmospheric [26] data, neutrinos are massive parti- 
cles. On the other hand, each of the extra Higgs scalar 
multiplets that we will introduce has already been con- 
sidered, separately, as viable extensions of the minimal 
SU(2) L ® U(1) Y model [27-29]. However, in general the 
doublets can be coupled to all fermion doublets, hav- 
ing flavor changing neutral currents if we will. In other 
words, there is not an underlying reason to constraint the 
interactions of the multi-Higgs extensions. Below we will 
show how this reason arises in an extension of the DFS 
axion model. 

The representation content is the following: Ql = 
(ud)l ~ (2,1/3), L L = {vl) T L ~ (2,-1) denote any 
quark and lepton doublet; ur <~ (1,4/3),(2r~ (1,-2/3), 
Ir <~ (1,-2), z/r ~ (1)0) are the right-handed compo- 
nents; and we will assume that each charge sector gain 
mass from a different scalar doublets, i.e., $„, <&d, 
and <!>,, generate Dirac masses for tt-like, d-like quarks, 
charged leptons and neutrinos, respectively [all of them 
of the form (2, +1) = (ip + , f°) T ]- We also add a neutral 
complex singlet </> <~ (1, 0) as in [7], a singly charged sin- 
glet h + <~ (l,+2), as in the Zee's model [27] and finally, 
a triplet T ~ (3, +2) that can be written as [28] 

where e = ir 2 . 

Next, we will impose the following Z13 symmetry 
among those fields: 

Q — > lo 5 Q, u R — > lo 3 u r , d R — > ui^dR, 
L — > uj e L, vr — > ujqVr, Ir — > uj 4 Ir, 
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(2) 

with w fe = e M / 13 , fc = 0, 1, 6. 

With this representation content and the Z13 symme- 
try defined in Eq. (2) the allowed Yukawa interactions 
are given by 

-C Y = Qi L (F ia u aR $ u + F ia d aR $ d ) +T^E(G ab u bR ^ v 
+ Gabl bR $t) + f ab TL~LrL bjL e l0 h+ + H.c, (3) 

where $ = e$*, with i = 1,2,3 and a = e, /z,r; G, G, F, F 
are arbitrary 3x3 matrices and f ab is a 3 x 3 antisymme- 
tric matrix. Notice that we have the Zee's interaction [27] 



but the interaction (Li a i y ) c F ab (eT ■ T)ijL b ^j is not al- 
lowed by the Z13 symmetry, i.e., this is not an automatic 
symmetry of the model. Moreover, the Z13 symmetry 
also implies that there is no flavor changing neutral cur- 
rents since there is a doublet for each charge sector in 
Eq. (3). The neutrino interactions with the doublet $„ 
imply a Dirac bare mass; while the singlet h + induces, 
through the radiative corrections, like in the Zee's model, 
Majorana masses to the left-handed neutrinos [30]. 

The most general scalar potential among the several 
scalar multiplets introduced above must be invariant un- 
der the gauge SU (2)l®U(1)y and also under the discrete 
Z13 symmetry. However, an extra discrete Z3 symmetry 
will be added in order to avoid an interaction in the scalar 
potential and in Eq. (3). We then have 



V = lil$>\$k + 4<f>*<f> + MrTr(TtT) + fih+h,- + \ k (<p{ <5> k ) 2 + \ kk ' \${$ k <\ 2 + A fe< ^$ fe )0* + \f^ k Tr(T^T) 
+ \f '($£Tt)(T* fc ) + \ kh *\$ k hrh + + \i(<f>*<f>) 2 + A 2 (TrTtT) 2 + \ 3 Tr(T^T) 2 + A 4 Tr(T t T)|0| 2 + A 5 (/T/i + ) 2 
+ X 6 T r (T^T)h-h+ + [Ai$J*,*l;$, + A 2 ^T$„<£ + A 3 $t $„$t $ ; + f&faf + / 2 $^T$„ + / 3 $f T$„ 
+ H.c], (4) 



where fc 7^ fc' run over all doublets, i.e., fc, fc' = u, d, I, v 
and we have omitted summation symbols. The trilincar 
§ u e§ih~(f>* is allowed by the Z13 symmetry but it is for- 
bidden by the Z 3 symmetry with parameters denoted by 
Wo, tDi, and Cj^ 1 if transform with Cu^ 1 and v R , l R 
transform with u>i while all the other fields remain in- 
variant i.e., transform with tDp. This Z3 symmetry for- 
bids also a Majorana mass term iy aR ) c M ab v bR in Eq. (3). 
It is in this context that the PQ symmetry is an auto- 
matic symmetry of the classical Lagrangian of the model 
and the axion mass and the 6 angle are protected against 
gravitational effects, as discussed in the Introduction. 

We have confirmed that the axion is an invisible one. 
After redefining the neutral fields as usual <I>° = (v k + 
Re$g + zIm$°)/V2, T° = (v T + RcT° + TmT°)/\/2, and 
4> = {vtf, + Kccj) + iIm(j))/y/2 there are only two neutral 
Goldstone bosons: one which is eaten by the Z°, G°, 
and another one which has a small mixture with the G°, 
the axion which is almost singlet, i.e., a ~ Im</>. This has 
been done considering the full 6x6 mass matrix in the 
pseudoscalar sector. 

Besides, the fact that the Peccei-Quinn chiral symme- 
try is an automatic symmetry of the model, the discrete 
Z13 symmetry suppresses also gravitational effects up to 
operators of dimension twelve. It means that the contri- 
butions to the axion mass and to the 9 angle are like those 
we have shown above. The assignment of the PQ charges 
for all the fermion and scalar fields in the model are the 
following (using the notation ip' — e""*^ ip with the 
PQ charge of the multiplet ip): Xq = —3, X UR = 3, 
Xd R = 3, Xl = —6, Xi R = —8, X VR = 4, X<s> u = 6, 



X^ d — — 6, X^ l — 2, X<$ u — 10, Xx — 12, Xj t + — 12, 
X<p = —12. It is straightforward to verify that the 
Yukawa interactions in Eq. (3) and the scalar potential 
in Eq. (4) are invariant under the PQ transformations 
given above. In the present model we have still a global 
discrete symmetry, namely Zi% C [/(I)pq. Hence, the 
model has the domain wall problem [31] and the usual 
mechanism for avoiding it can be applied [32,33]. 

Larger Zn symmetries are possible if we add more 
fields in any electroweak model. However, notice that 
if N is a prime number the axion can transform under 
this symmetry with any assignment (but the trivial one); 
otherwise, we have to be careful with the way we choose 
the singlet <ft transforms under the Zn symmetry. 

In the present model the axion couples to quarks and 
leptons (including neutrinos). Those couplings are pro- 
portional to X fmfv^ 1 , where Xj and m/ are the PQ 
charge and the mass of the fermion /. On the other 
hand, the axion-photon coupling c a77 is, in general, pro- 
portional to 2iVj^y XfQ"j — 1.95, where the sum is 
over all generations and the numerical factor comes from 
the light quarks through a Fujikawa relation [34] and it 
is the same in the present model; -/Vdw is the domain 
wall factor which is equal to 18 in the present model. We 
have Ca-yj ~ —0.62, and not 0.75 as in the invisible axion 
of Rcfs. [7,35,36]. The coupling of the axion with neutri- 
nos in the present model may have astrophysics and/or 
cosmological consequences. 

It is still possible to assign the PQ charge as in Sec. II 
of Rcf. [22] and the PQ charges are quantized. However, 
in that case it appears as an extra Goldstone boson cou- 



3 



pling with leptons and photons. Thus, in order to avoid 
such an extra Goldstone we have to break softly the Z3 
symmetry, say, by adding /z 2 $;$„ in the scalar potential. 

Summarizing, we have obtained a model with an au- 
tomatic PQ symmetry, with an invisible axion which is 
protected against gravitational effects in the context of 
a multi-Higgs extension of the standard model. This is 
important for three reasons. First, multi-Higgs exten- 
sions, with or without supcrsymmctry, are used by the 
experimentalists as reference models for searching extra 
scalar particles [29,37]. Second, these models are used 
by phenomenologists to get an appropriate mass matrix 
in the lepton and quark sectors. In fact, as we said be- 
fore, all the extra fields that we have added to the min- 
imal version of the model have already been introduced 
in the literature with different goals; we have just put 
all of them together. The necessity of protecting the in- 
visible axion from gravitational effects, as we have done 
above, gives a rationale for those scalar multiplcts. Fi- 
nally, we would like to mention that one of the issues 
that remains arbitrary in most multi-Higgs extensions is 
the flavor violating neutral currents in the scalar sector. 
We can implement it or not at will. However, in our case 
these effects are automatically avoided since the Yukawa 
interactions are already determined for the Z13 symme- 
try. Then, our model automatically does not have flavor 
changing neutral currents in all sectors of the model. 
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